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i Abstract. INTEGRAL regularly scans the Galactic plane to search for new objects and in particular for absorbed sources with 

. the bulk of their emission above 10 - 20 keV. The first new INTEGRAL source was discovered on 2003 January 29, 0.5° from 

' the Galactic plane and was further observed in the X-rays with XMM-Newton. This source, IGRJ1 63 18-4848, is intrinsically 

strongly absorbed by cold matter and displays exceptionally strong fluorescence emission lines. The likely infrared/optical 
counterpart indicates that IGRJ1 63 18-4848 is probably a High Mass X-Ray Binary neutron star or black hole enshrouded in a 
Compton thick environment. Strongly absorbed sources, not detected in previous surveys, could contribute significantly to the 
Galactic hard X-ray background between 10 and 200 keV. 
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1. Introduction of 508 ksec. The position of the source (RA = 16 h 31.8 m 

• ' / , , i . . and DEC = -48°48') was determined with an accuracy of 

k^l X-ray binaries (where the compact obi ect is a neutron star or or Trn TU . 1C aoao • a + ±a ± on i \t vu 

^< , - / , , 11 . i 2 • IGR J16318-4848 is detected up to 80 keV with a mean 

r . ^ black hole) can become strong hard X-ray s emitters when ac- ork crk , w fl f . / in _n -2 -l c- + ^ c \ 

./i i . . 20 - 50 keV flux of 6 x 10 11 erg cm z s . Significant (> 5<x) 

& cretion takes place. Among the -300 known X-ray binaries in . , . . . . . T u . , nrkrk 

. - r b . J intensity variations occur on time scales as short as 1000 sec. 

our Galaxy and the Magellanic clouds, a few systems show 

strong intrinsic photo-electric absorption: GX 301-2 (Swank We selected 70 pointings of the INTEGRAL core program 
et al. 1976), Vela X-l (Haberl & White 1990), CI Cam (Boirin between revolution 36 and 50 in which IGR J16318-4848 was 
et al. 2002). Moderate absorption was also detected in a few located at less than 5° from the center of the field of view. 
X-ray bursters (Natalucci et al, 2000). We report here on the since IGR J16318-4848 is often below the detection level, fur- 
discovery of IGR Jl 63 18-4848, a Compton thick X-ray binary ther selection was applied on the luminosity of the source, 
in which the X-ray obscuring matter has a column density as We selected pointings where the source was detected by the 
large as the inverse of the Thomson cross section. IBIS analysis software (Goldwurm et al. 2003) and out of 

them only 8 pointings (from 2003 March 3 to 2003 March 
14) for which we were able to produce spectra of sufficient 

2. High energy Observations and data analysis quality. This resulted in a 17.5 ksec average spectrum that 

T/-TD ti^io AO AO a- a • +u tatt^d a t • is displayed in Fig. El The current ISGRI energy correction 

IGR J163 18-4848 was discovered using the INTEGRAL lm- , . . ^ , 

tt> jc/jc^ r> j atu -4. i onno T u * i onm\ and response matrices are preliminary. Throughout this anal- 
ager IBIS/ISGRI (Ubertini et al. 2003; Lebrun et al. 2003) on . r Tr ^„ T h i i 

T on \ i i i ysis we used an ISGRI ancillary response that was modified 

2003 January 29 (Courvoisier et al. 2003a) and was regularly J , . , _ , J ^ r t _ , ^ T , , _ 

i i n . . i -|— 1« m i 1C Ari t xr to obtain a good fit to the spectrum or the Crab Nebula. The 

observed for two months. Figure 111 shows the 15 - 40 keV _ & , , , _ n A n A n 

TC „ DT , . A4.u -4u i ISGRI spectrum of IGR J16318-4848 could be represented by 
ISGRI sky image around the source with an accumulation time , . t t . t ^ ^ „+i ? i „ 
a power law with a photon index, Y = 2.7 and a flux 

Send offprint requests to: Roland.Walter@obs.unige.ch F20-100 keV = 1.6 X 10" 10 erg cm" 2 sec" 1 . 
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Fig. 1. INTEGRAL ISGRI 15-40 keV sky image of the Norma region in Galactic coordinates. This image accumulates 508 ksec 
of INTEGRAL core program data. The red circle indicates IGR J16318-4848. 



IGRJ16318-4848 was observed for 28 ksec by XMM- 
Newton on 2003 February 10. A single X-ray source was found 
within the INTEGRAL uncertainty circle in the EPIC PN and 
MOS cameras (Striider et al. 2001; Turner et al. 2001) at a po- 
sition of RA = 16 h 31 m 48.6 s and DEC = -48°49 , 00" with 
a 4" uncertainty (Schartel et al. 2003). The X-ray spectrum 
was immediately recognized as exceptional featuring strong 
photo-electric absorption, the associated Fe absorption edge at 
7.1 keV and fluorescence line emission of mostly neutral Fe Ka 
(6.4 keV), Fe KJ3 (7.1 keV) and Ni Ka (7.5 keV). The XMM- 
Newton spectrum of IGR J16318-4848, presented by Matt & 
Guainazzi (2003) is significantly flatter than the ISGRI spec- 
trum. 

All spectral uncertainties are given at 90% confidence 
for a single interesting parameter, unless indicated otherwise. 
The abundances of Anders & Grevesse (1989) and the photo- 
electric cross section of Verner et al (1996) are used throughout. 

To analyze the ISGRI and EPIC data together we extracted 
EPIC spectra using version 5.4. 1 of the XMM Science Analysis 
System (SAS) software. The data were first screened for en- 
hanced variable background by filtering out the time intervals 
where the count rate above 10 keV was higher than the thresh- 
old count rate (18 for MOS2 and 60 for PN per 100 sec bin). 
The total exposure after screening resulted in 24 ksec for MOS2 
and 21 ksec for PN. Source events were subsequently extracted 
from a 25 " radius circle centered on the source. A second circle 
with the same radius was fixed on a comparable region on the 
detector to serve as background. Standard SAS tools were used 
to calculate the instrumental response and the effective area for 
the extracted spectra. 



Table 1. Best fit parameters (90% confidence). All parameters 
are free in model 1. Cisgri is fixed to 0.67 in model 2 



Parameter 


Model 1 


Model 2 


Unit 




X 2 /d.o.f. 


333 / 322 


336 / 323 






ClSGRI 


43+ - 14 

-0.23 


0.67 (fixed) 








1.96 + 0.07 


2.07 + 0.10 


10 24 


cm 2 


r 


1.6 + 0.3 


1.97 + 0.17 






Pwl h keV 


0.08 + 0.03 


0.12 + 0.03 


ph/k 


ev cm 2 s 


Fe abs 


0.82 + 0.05 


0.79 + 0.04 


Z Q 




Fe Ka E. 


6.405+0.003 


6.405 + 0.003 


keV 




Fe Ka Flux 


1.84 + 0.09 


1.75+0.06 


io- 4 


ph/cm 2 s 


Fe KJ3 E. 


7.07 + 0.01 


7.07 + 0.01 


keV 




Fe K/3 Flux 


0.32 + 0.08 


0.30 + 0.07 


io- 4 


ph/cm 2 s 


Ni Ka E. 


7.46 + 0.02 


7.46 + 0.02 


keV 




Ni Ka Flux 


0.09 + 0.03 


0.09 + 0.03 


io- 4 


ph/cm 2 s 


Intr. line width 


<20 


<20 


eV 





The simultaneous fit of the EPIC PN and MOS2 and ISGRI 
spectra of IGR J163 18-4848 using an absorbed power-law 
continuum, free Fe abundance and three Gaussian emission line 
model (model 1) resulted in a reduced^ 2 of 1.04 for 322 de- 
grees of freedom (d.o.f.). A normalization constant Cisgri was 
introduced in the model and was left as a free parameter in the 
fit in order to account for cross calibration uncertainties and 
the non simultaneity of the ISGRI and EPIC observations. The 
best fit parameters, listed in table [I] are compatible with those 
of Matt & Guainazzi (2003). The powerlaw flux in table 1 is 
corrected for the effects of photoelectric absorption but not for 
the effects of Compton scattering. The line fluxes are observed 
fluxes, not corrected for any absorption effect. 
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Fig. 2. The EPIC PN, EPIC MOS2 and ISGRI photon spectra of IGR J16318-4848 along with the best fit model and residuals 
for Cisgri fixed to 0.67. The fit gave a reduced^ 2 of 1.04 for 323 d.o.f. (table 1). 



Examination of the residuals of model 1 shows that the 
spectral model is systematically flatter than the ISGRI spec- 
trum. In addition as the ISGRI spectrum represents the high 
state of the source and the EPIC data correspond to a mix of 
different levels we would expect Cisgri to be > 0.67, the value 
derived using our ISGRI response matrix and the EPIC/ISGRI 
inter-calibration obtained on 3C273 (Courvoisier et al. 2003b), 
in contrast to the best fit value of 0.43. 

Matt & Guainazzi (2003) noted that the spectral slope was 
weakly constrained by the EPIC data alone as it correlates with 
the absorbing column density. Cisgri correlates with the power- 
law index as well. Constraining Cisgri to be > 0.67 gives T = 
1.9 ± 0.2 and a column density N n = (2.1 ± 0.1) x 10 24 cm -2 . 
The photon index is consistent with that derived from the fit to 
the ISGRI spectrum alone. Table [I] gives the best fit parameter 
for a model with Cisgri fixed to 0.67 (model 2), the photon 
spectrum and the residuals are shown in Fig.|2| 

An exponential cutoff power-law continuum model also fits 
the data well providing a x 2 of 328 for 322 d.o.f. with a flat 
power-law and a cutoff energy of 15 ± 5 keV. Simultaneous 
spectral observations performed below and above 20 keV are 
necessary to constraint better the spectral model. 

The observed ratio of the intensities of Fe Kfi and Fe Ka 
(0.17 ± 0.09) is consistent with the expected value of 0.14 
(Kaastra & Mewe 1993). The intensity ratio between Ni Ka 
and Fe Ka of 0.05 ± 0.02 is close to the solar abundance ra- 
tio of 0.03-0.045 (Molendi et al, 2003). The strength of the Fe 



edge at 7.1 keV also corresponds to that expected from a solar 
abundance. 

The centroids of the Fe Ka and KJ3 lines corresponds to Fe 
that is ionized between 2 to 6 times (l<x level) (House, 1969) 
and to a ionisation parameter S « 0.05 erg cm s" 1 that is ex- 
pected to be variable across the absorbing matter. The position 
of the Fe absorption edge also indicates the presence of Fe that 
is ionized less than 2 times. Note that the systematic uncertainty 
on the line and edge energies is 10 eV. 

In contrast to Matt & Guainazzi (2003), we did not find that 
the presence of a Compton shoulder to the Fe Ka line was re- 
quired by the data. This could be related to data selection and 
reduction. A firm detection of the Compton shoulder should be 
confirmed by future observations. The absence of detection of 
the Compton shoulder is however consistent with the conclu- 
sion of Matt & Guainazzi that the average Nu on which fluo- 
rescence takes place could be smaller than that on the line of 
sight. We consider the flux of the Compton shoulder derived by 
Matt and Guainazzi (2003) as an upper limit. 

We also did not find evidence for an excess of emission 
below 5 keV which could be explained by the steeper contin- 
uum derived from the INTEGRAL data when compared with 
the use of XMM data alone. Our data do not therefore show ev- 
idence for reflection as could be inferred from the low energy 
excess. A pure transmission geometry, with in-homogeneously 
distributed absorbing mater is sufficient. 
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Matt & Guainazzi (2003) noted that the source flux vari- 
ations were intrinsic to the source (and not related to absorp- 
tion). We extracted X-ray light curves for the continuum and 
for the Fe Ka line (Fig.0) using time bins of 200 sec. During 
the fast continuum rise (around time 19000 sec in Fig.0 the 
count rate varies significantly in 200 seconds and a cross corre- 
lation analysis between the 6.2 - 6.6 keV light curve (strongly 
dominated by the Fe Ka line) and the 8 - 12keV continuum 
light curve indicates that any differences in mean arrival time 
of the continuum and Fe Ka line emission are smaller than 200 
sec. The ratio between the Fe Ka and continuum shows how- 
ever some significant variations which, as pointed out by Matt 
& Guainazzi (2003), could indicate variations of the properties 
of the cold matter on time scales of 10 4 sec. 
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Fig. 3. EPIC PN light curves accumulated in different energy 
ranges for IGR J 163 18-4848. The Fe Ka light curve has been 
accumulated between 6.2 and 6.6 keV. The 8-12 keV light 
curve is unaffected by line emission. 



3. Discussion 

3.1. X-ray variability 

The variability time scale and the maximum delay observed 
between the Fe Ka line and the continuum variations limit the 
size of the zone in which fluorescent emission takes place and 
its distance to the X-ray source to 10 13 cm. It is therefore very 
unlikely that IGR Jl 63 18-4848 is an extragalactic source such 
as a Seyfert II galaxy or an Ultra Luminous Infra Red Galaxy as 
the width of a fluorescence line (given by the Keplerian veloc- 
ity) emitted at a distance of 10 13 cm would be orders of magni- 
tude larger than observed. 

Following the INTEGRAL discovery of IGR J16318-4848, 
a re-analyis of archival data showed that IGR J 163 1 8-4848 had 
been weakly detected in 1994 September by ASCA (Murakami 
et al. 2003), with similar flux and Nu as observed in 2003. The 
ASCA spectrum suggested the presence of a strong Fe Ka line 
(Revnivtsev et al. 2003). However, the Beppo-SAX Wide Field 
Camera (WFC), that observed the field almost continuously for 
6 months each year between October 1996 and 2002, has never 
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Fig. 4. EPIC MOS, 2MASS (J, H, K) and DSS images of the 
counterpart of IGR J163 18-4848. All images have the same 
scale. The XMM error circle (4" radius) is shown on each im- 
age. 



detected the source. This indicates that it was -10 times fainter 
on average during those periods (In 't Zand 2003). This suggests 
that IGR J16318-4848 was active in 1994 and 2003 and that it 
has been quiet over periods of many months. It is remarkable 
that the flux and the absorption observed in the active states in 
1994 and 2003 are very similar. 

3.2. Counterpart 

A possible counterpart to IGR J16318-4848 in the Digitized 
Sky Survey (DSS-II/USNO-B1.0), Two Microns All Sky 
Survey (2MASS), and Midcourse Space Experiment (MSX) 
data was reported by Foschini et al. (2003). The respective im- 
ages are shown in Fig. @] In the 2MASS, the counterpart has 
/ = 10.2, H = 8.6, K = 7.6 and an uncertainty of ±0.3 mag. It 
is also clearly detected in the I band of the second DSS and also 
in the R band (USNO-B 1). The R magnitude is reported to vary 
over an interval of 50 years between 17.3 ± 0.3 and 18.4 ± 0.3. 
The flux density in the MSX A band is 0.46 Jy. 

Radio observations were performed with the Australia 
Telescope Compact Array (ATCA). Observations have been 
conducted at 4.8 and 8.6 Ghz (with a total bandwidth of 128 
MHz) on 2003, February 9, starting at 18:00 and finishing at 
04:00 the following day, with a total of 1.33 hours (spread 
over the full ATCA run) on IGR J16318-4848. No source was 
detected with a 1 cr upper limit of 0.1 mJy both at 4.8 and 
8.6 GHz. 

The near infrared spectral energy distribution corrected for 
the effect of various possible values of the galactic absorption 
(Lutz et al. 1996) was analyzed. The visual extinction factor 
derived from the galactic TVh along the line of sight is Ay = 11. 
The absorption could however be smaller if the source is nearby 
or larger as the radio measurements have limited spatial reso- 
lution. The maximum reddening compatible with an infrared 
spectrum not flatter than a black body spectral component 
(Ay = 20) is two orders of magnitude smaller than the absorp- 
tion observed in the X-rays. We conclude that the source of 
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the IR/optical emission is located at > 10 13 cm from the X-ray 
source. 

If the reddenning is low (Ay < 11) the infrared compan- 
ion of IGR J 163 18-4848 is a low mass red giant star with a 
luminosity of 100 d^ c L©. If the reddening is strong (Ay « 
20) the companion is a massive supergiant star of luminosity 
10 5 <^ kpc L©. In the rest of this discussion we will assume that 
the infrared/optical counterpart and IGR J 163 1 8-4848 are a bi- 
nary system. This however remains to be verified. 

If the system is a Low Mass X-Ray Binary (LMXRB) lo- 
cated at 1 kpc, the unabsorbed 20 - 50 keV luminosity of 3 x 
10 34 erg s" 1 corresponds to the quiescent state of those systems. 
In that state, neutron star systems emit thermal emission (kT « 
100 eV) that dominates the hard tail below 3 keV (Rutledge 
et al. 2002) which is not detected in IGR J16318-4848. Black 
hole LMXRB systems in quiescence have lower luminosities 
(Kong et al. 2002) and do not show strong absorption. 

Alternatively, if IGR J16318-4848 is a High Mass X-ray 
Binary (HMXRB) located at 5 kpc, the unabsorbed 20 - 50 keV 
luminosity of 7x 10 35 erg s" 1 indicates that moderate (wind) ac- 
cretion is taking place. Strong Fe Ka lines have been observed 
in other HMXRB systems. In Vela X-l, large line equivalent 
widths were observed during eclipses when only scattered X- 
rays are observed, the absorbing Nu remained however smaller 
than observed in IGR J16318-4848 (Pan et al. 1994). In GX 
301-2, the Nu and the equivalent width of the Fe Ka line are 
variable and correlated. White & Swank (1984) explained this 
behavior by variations of the stellar wind velocity and den- 
sity along the orbit of the compact source. The Fe Ka equiva- 
lent width and the N n observed from IGR J16318-4848 match 
the extreme of the correlation found in GX 301-2. This sug- 
gests that the stellar wind accreting onto the compact source of 
IGR J163 18-4848 could form a dense spherical shell in which 
fluorescence and absorption takes place. This possibility was 
also proposed by Revnivtsev et al. (2003). Note that the expo- 
nential cutoff powerlaw model that was used to represent the 
data is rather typical for accreting X-ray pulsars. 

The ionisation parameter derived from our observations can 
be used to estimate the distance between the fluorescing mate- 
rial and the X-ray source D = (a/D) L/E Nu = 10 13 (a/D) cm 
where a is the shell thickness. This distance is also compatible 
with the line variability. This distance compares better with the 
companion star radius than with the accretion radius (unless the 
stellar wind velocity is small). 

We searched unsuccesfully for the presence of pulsations 
that would be a clear signature of a neutron star in the system. 
The counting rate of the source during the XMM observation 
does not allow any definitive answer as the upper limit on the 
relative amplitude for a rather broad 0.15 Hz quasi periodic os- 
cillation (QPO) (e.g. Vcentroid/FWHM = 3) is ~ 30% at the 
3 cr confidence level, which is not constraining. An 8000 sec 
observation with the Proportional Counter Array (PCA) on- 
board RXTE did not find evidence for QPO either (Swank & 
Markwardt, 2003). 



3.3. Contribution to the X-ray background 

Absorbed sources could contribute significantly to the Galactic 
high energy background. Valinia et al. (2000a) modeled the 
Galactic background emission observed by the PCA on RXTE 
and OSSE on CGRO using a specific spectral component dom- 
inating between 10 and 200 keV. The flux of that component 
(2.5 x 10" 11 erg cm" 2 s" 1 deg" 2 at 40 keV) is variable indicating 
that an important fraction of that emission comes from point 
sources. The spectral shape of IGR J163 18-4848 roughly cor- 
respond to the empirical spectral model used to represent the 
Galactic background. The intensity of the background emis- 
sion could be explained by 0.4 sources per square degree at the 
level of IGR Jl 63 18-4848 in a band of a width of few degrees 
around the Galactic plane. 

The hard X-ray background also displays a strong Fe Ka 
emission line which is attributed to the thermal diffuse emission 
dominating below 10 keV (Valinia et al. 2000b). A population 
of intrisically absorbed sources similar to IGR J163 18-4848 
would only contribute at a level of 10% to the background line 
intensity. 

In spite of the BeppoSAX WFC long term monitoring of 
the Galactic center only few strongly absorbed sources have 
been discovered so far (Ubertini et al. 1999). INTEGRAL is able 
to make sensitive search for highly absorbed sources (Lebrun 
et al. 1999) such as IGR J16318-4848, IGR J16320-4751 
(Rodriguez et al. 2003) and IGR 16358-4726 (Revnivtsnev et 
al. 2003). Figure Q] shows the detection by INTEGRAL of one 
point source per degree of Galactic longitude in the Norma re- 
gion. Those sources could indeed explain a significant fraction 
of the Galactic diffuse emission. 
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